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Abstract 



• A random phase property is proposed for products of random matrices drawn from 

one of the classical groups associated to the 10 Cartan symmetry classes. It allows to 
calculate the Lyapunov spectrum explicitly in a perturbative regime. These results apply 
to quasi-one-dimensional random Dirac operators which can be constructed for each of the 
symmetry classes. For classes corresponding to quantum Hall systems, quantum spin Hall 
systems and Z2-topological superconductors the random Dirac operators have vanishing 
" Lyapunov exponents and almost surely an absolutely continuous spectrum. 

m- PACS: 72.15.Rn, 73.23.-b 

o 

CN ; The theory of products of random matrices has successfully been applied to the analysis of 

localization effects in quasi-one-dimensional disordered systems. In the mathematical literature 
[BLj . particular focus has traditionally been put on time-reversal invariant systems without 
spin, for which the transfer matrices are in the symplectic group. On the other hand, it is well 
known from work of Altland and Zirnbauer |AZ] (see also [SRFLl iRSFLj ) that there are also 
disordered systems of non-interacting Fermions with time-reversal, sublattice (also called chiral) 
and/or particle- hole symmetries (denoted by TRS, SLS and PHS below). In dimension 1 this 
leads to transfer matrices in other classical groups (see Section [T] below for terminology), and 
actually each of the non-exceptional groups of the Cartan classification appears as the transfer 
matrix group of an adequate quasi-one-dimensional model. For example, an SLS appears for 
lattice operators typically at special energies when there is a sublattice structure. The best 
known example is probably the Dyson model with off-diagonal disorder, for which the transfer 
matrices at zero energy have such a chiral symmetry (it is in the chiral orthogonal class). 
On the other hand, the Bogoliubov-deGennes (BdG) classes describe quasi-particles in dirty 
superconductors. These quasi-particles can be either particles or holes, and the Hamiltonians 
then have a supplementary PHS. Random Dirac operators with adequate physical symmetries 
allow to construct concrete models having transfer matrices (fundamental solutions) in each 
class |BMSAt IBFGMj . The interplay between physical symmetries of these Hamiltonians and 
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the symmetries of the transfer matrix group are reviewed in Sections [3] and [4] below. Of course, 
apart from this physically motivated application the study of random products in each of the 
classical groups is of intrinsic mathematical interest as well. 

The most interesting quantities associated to such random products are the Lyapunov ex- 
ponents, which are the exponential growth rates of their singular values. The collection of 
these exponents constitute the Lyapunov spectrum. In the physics literature, the approach 
most widely used for their study is based on the Dorokhov-Mello-Pereyra-Kumar (DMPK) 
equations [Port IMPKj . These equations are universal (in the sense that they do not contain 
model dependent information other than the above mentioned symmetry properties), but their 
derivation relies on the so-called maximal entropy Ansatz. A recently proposed alternative 
approach [RS2j is closer to the theory of products of random matrices and the numerical trans- 
fer matrix method. One considers the random action of the transfer matrices on the compact 
Grassmannian. This leads to a Markov process and then one supposes that the random phase 
property (RPP) holds, namely that the invariant measure of the Markov process is given by 
geometric invariant measure on the Grassmannian. In [RS2j the RPP was only studied for the 
three standard symmetry classes, and here this notion is extended to treat the other symme- 
try classes of BdG and chiral type. Models for which the RPP holds exactly can easily be 
constructed, see Section [2 As discussed in |RS2j . the RPP is considerably weaker than the 
maximal entropy Ansatz and it is susceptible to numerical verification in concrete models. 

Whenever the RPP holds, the calculation of the Lyapunov spectrum reduces to the eval- 
uation of the averages of certain functions over the invariant measures on the Grassmannian. 
These invariant measures are given in terms of the Haar measure on the maximal compact 
subgroups of the given classical group, and all these compact subgroups are in turn given in 
terms of the unitary, orthogonal or symplectic group. While in principle it is known how to 
calculate these averages [Weil ICS} ICStj . it becomes feasible only in a perturbative regime of 
small randomness. Then one only needs the second and some fourth moments of these Haar 
measures. The necessary formulas are collected in the Appendix. The outcome of the rather 
tedious calculations are nice and compact formulas for the Lyapunov spectrum in a perturba- 
tive regime. They are listed for each symmetry class separately in the subsections of Section [2J 
The Lyapunov exponents 7 P > 7 p+ i in a system with N channels and a centered random 
perturbation of strength are in all cases roughly given by the same formula 

X2 N-7] P 3 
7 P A h U ( X ) ' 

with rj = 1 or n = 2, but particular focus is on calculating the prefactors in this formula in each 
symmetry class. This also allows to calculate the smallest non-negative Lyapunov exponent 
which in physical models is interpreted as the inverse localization length. The precise form 
of the prefactors is particularly important for the study of systems with a small number of 
channels. These general results on the Lyapunov spectrum apply directly to the model Dirac 
operators of Section HJ 

It is possible that the products of random matrices have vanishing Lyapunov exponents 
which are symmetry enforced, namely they result from the defining relations of the classical 
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group considered. The most basic example of this type are products in the group U(iV, M) 
with M > N which have at least M — N vanishing Lyapunov exponents. This corresponds to 
a quantum Hall system with N — M edge channels (the quasi-one-dimensional operator then 
gives an effective description of the edge states). Subclasses correspond to the thermal and 
spin quantum Hall effect. Other interesting cases concern quasi-one-dimensional Hamiltonians 
in the Cartan class All and Dili which describe boundary states in two-dimensional quantum 
spin Hall systems and time reversal invariant dirty superconductors respectively. The parity of 
their channel number N can be understood as a Z 2 -invariant and, if it is non-trivial (namely, 
N is odd), then there are again at least two vanishing Lyapunov exponents. This leads to 
derealization and to absolutely continuous spectral measures. The case All was studied in 
detail in [SS2j . Here in Section [5] we discuss how these results transpose to the case Dili as 
well as quantum Hall systems. 

1 A list of the classical groups 

The classical groups in the sense of Weyl incorporate certain families of matrix groups. The 
list is given in the first chapter of |GW] . There is also a close connection to the Cartan 
classification of symmetric spaces because each classical group from the list below corresponds 
to a non-compact Lie algebra in duality with a classical compact symmetric space |Hel] . As 
such there is a Cartan label associated to each classical group and this label has been used 
heavily in the physics community \AZ\ IBFGMj . For this reason, we also choose to attach the 
Cartan label to each classical group Q, but at the same time also the classical terminology is 
recalled in each case below. Because it will be heavily used later one, we also recall the maximal 
compact subgroup U for each classical group Q as well as peculiarities of the spectral theory 
in Q. The groups Q are listed in Table 1, albeit in a different order than we go through them 
here. The groups will be described using the matrices 




where the identity 1 is of adequate size, mostly an iV x N matrix. These are real matrices 
satisfying J* = J, K* = K, I* = —I and 1 = J 2 = —I 2 = K 2 which is why we call J and 
K an even, but I an odd symmetry. They do not commute, and one has K J = I. Of course, 
one could use Pauli matrices instead. Writing i = y/— 1, the three involutions J, K and 1 1 are 
isospectral, and one can pass from one to another by a unitary transformation. In particular, 
the Cayley transformation C achieves the following 

C*JC = il , CJC* = K , C = 4= 

Moreover, C*IC = iK, C l KC = 1, CKC 1 = iJ and CJC* = -il. Among the possible 
realizations of the classical groups we choose ones for which the maximal compact subgroup 
takes a particularly simple form, since this will allow to compute averages over these subgroups 
later on. 
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1.1 Class A 



In Class A there are no further symmetries imposed. Thus Q A = GL(A r , C) is the set of all 
invertible N x N matrices. The maximal compact subgroup U A = U(iV) is given by the unitary 
matrices. 

1.2 Class AI 

Class AI contains the subgroup of Q A composed of the real matrices (invariant under the 
operation of complex conjugation): 

g A1 = [T e GL(iV,C) |T = T } = GL(iV,R). 

Clearly the maximal compact subgroup U A1 = O(N) is composed of the orthogonal matrices. 

1.3 Class All 

The definition of Class All is 

gAix = | T e GL (2iV,C) | PTI = T } = U*(2iV) . 

Writing out the symmetry one readily checks that T G U*(2iV) is of the form T = (^) with 
entries given by iV x N complex matrices. One obtains an N x N quaternion matrix from 
T by setting T(T) = A + Bj where j is the second of the standard unit quaternions 2, j, k 
satisfying i 2 = j 2 = k 2 = ijk = —1. This map J 7 is a *-isomorphism from U*(2N) to the group 
GL(iV, M) of invertible quaternion matrices. Thus the classes A,AI,AII simply correspond to 
the invertible matrices over the field C, M, H respectively. The maximal compact subgroup of 
U*(2iV) is 

A*A + B t B = l, A*B t = B t l\ = SP(2N) . 

The spectral theory of T G U*(2iV) has the following feature. If Tv = Xv for some v G <C 2N and 
A G C, then T Iv — XIv and furthermore the vectors v and Iv are linearly independent. Indeed, 
suppose v = m satisfies v = filv for some non- vanishing [i G C. Then a = —jjb and b = fia, 
which combined imply a = — \fi\ 2 a so that a = b = 0. This implies that the real spectrum of 
T G U*(2N) has even geometric multiplicity, a fact that is called Kramers' degeneracy in the 
physics literature. Actually, all the above also holds for all (possibly non-invertible) matrices 
in the Lie algebra u*(2iV) (which by the above map J 7 is mapped to all N x N quaternion 
matrices) . 

1.4 Class AIII 

The group of Class AIII is defined by: 

gAiu = | T 6 GL ( 2 jV, C) \ T* JT = J} = U(iV, N) , 



(2A0 n U(2iV) = j fj 
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for which maximally compact subgroup is 



V, W e U(N)j = U(N) x U(JV) . 

Elements of U(N, N) are often also called J-unitaries. Their spectrum has a reflection property 
around the unit circle S 1 because T — XI = AJ(T*) -1 (T — (A) _1 1)* J so that the invertibility of 
T — XI is indeed equivalent to the invertibility of T — (A) -1 1. Note, however, that in general the 
two eigenvectors satisfy no relation so that one also cannot conclude that each eigenvalue on the 
unit circle is degenerate. Another important feature of the spectral theory in U(iV, N) is that the 
eigenvectors v and w to eigenvalues A and \i with /ZA 7^ 1 are J-orthogonal, namely w* Jv = 0. 
Indeed, Tv = Xv and Tw = fiw imply w*Jv = {n~ x Tw)* JX~ x Tv = (jlX)~ 1 w*Jv which implies 
the claim. In particular, this shows that eigenvectors v corresponding to eigenvalues off the 
unit circle are J-isotropic, that is satisfy v* Jv = 0. 

Class AIII also contains a generalization of the group Q Ain defined above: 

U(JV,M) = {T e GL(N + M, C) \ T* GT = G } , G = I 1 " ° J , (3) 

\ u — ±m J 

where, say, M > N are possibly different. Then XJ(N, M) n U(JV + M) = U(JV) x U(M) is 
the maximal compact subgroup. All the above arguments about the spectral theory still apply 
to T e U(A r , M). In particular, the spectrum always comes in pairs A, (A) -1 and eigenvectors 
to eigenvalues A and \x with JlX 7^ 1 are J-orthogonal. This implies that there are at most 
TV eigenvalues off the unit circle, or, alternatively stated, at least M — N eigenvalues on the 
unit circle (here both iV and M — N give the algebraic multiplicity). Indeed, let v±, . . . , v p be a 
maximal set of linear independent generalized eigenvectors of eigenvalues of modulus larger than 
1. By the above, they then span a p-dimensional G-isotropic subspace, namely a subspace on 
which G seen as quadratic form vanishes. But the maximal dimension of G-isotropic subspaces 
is N, so that p < N. 



U AUI = U(N, N) n U(2iV) 



V 

w 



1.5 Class CI 

Class CI consists of the following subgroup of Q A ™ = U(N, N): 

g cl = [T e GL(2AT, C) I T*JT = J, K*TK = T}. 
Its maximally compact subgroup is 

u cl = g cl n u(2iv) = I K v e u(jv) j = u(jv) . 

Th group Q 01 is actually a disguised version of the real symplectic group, namely using the 
Cayley transform one has 

C*g CI C = {Te GL(2A^,C) j T*IT = I, T = T} = SP(2iV,M). 
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In this representation the maximal compact subgroup does not take such a simple form. The 
spectrum of T G Q CI has, of course, the S 1 -reflection property holding for any matrix in Q AIU , 
but, moreover, Tv = Xv leads to T Kv = XKv. Thus the eigenvalues come in quadruples 
A, A, A -1 , (A)" 1 . Let us point out that v and Kv are, in general, not linearly independent. 
Therefore one cannot conclude from the above that real eigenvalues are degenerate. 



1.6 Class Dili 

Also Class Dili is a subgroup of G AIU = \J(N, N): 

gum = | TG GL(2iV,C) \ T*JT = J, PTI = T} . 

Thus Q DIU = U(N, N) fl U*(2iV) and one can also see Q mu as a subgroup of Q AU . In particular, 
the map J 7 described in Section 11.31 is well-defined on Q DIU and the image J r (^ DIH ) is given by 
those quaternion matrices J-*(T) satisfying J r (T)*i jF(T) — 1 1. Its maximally compact subgroup 
is 

'V 0' 
V 



u mu = g mu n u(2iv) 



V e XJ(N) } ^ XJ(N) . 



Taking the Cayley transform of Q mu one finds, using C l IC — il, one of the classical groups: 

C*g mu C = {Tg GL(2iV,C) \T*IT = I, FTI = T} = SO*(2N) . 

Here the equation I*TI = T can equivalently be replaced by T l T = 1, because the relation 
T*IT = I holds. As to the spectral theory, again the S 1 -reflection property holds. Further- 
more, Tv = Xv implies T Iv — X Tv so that again the spectrum comes in quadruples as in 
qc\ ^ SP(2iV, K). Furthermore, because Q mn C Q An , the real spectrum has again a Kramers' 
degeneracy (see Section [T73|) . 

It is worth mentioning that the two symmetries J and I used in the definition of Q mi1 do 
not commute. If the matrix size is AN instead of 2N which we realize by tensorizing with C 2 
(or equivalently, iV is even in the above), then one can also use two commuting symmetries: 

(CByg DIII CB = {T G GL(4iV,C) |T*J®1T = J(8)1, 1 <g> I* Tl ® I = T} , 

where B = ^^J). Indeed, then B* = B~ l = B as well as B l l ® IB = -% I ® 1 and 
B*I ® IB = —I <S> 1. This representation of the group Q mu appears as transfer matrices in 
physical models with odd spin, where then 1 <S> I acts on the spin degree of freedom only (e.g. 
[RS2]). 

1.7 Class BDI 

Class BDI is yet another subgroup of Q AIU = XJ(N, N): 

g Bm = | Tg GL (2iV,C) \T* JT = J, T = T} = 0(N,N). 
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Of course, Q BUl is also the subgroup of Q AI = GL(2iV, R) specified by the relation T l JT = J. 
Its maximally compact subgroup is 



W m = g Bm n u(2iv) = { (o w) v > w G °( iV )} 



0(N) x O(N) . 



As matrices in 0(N,N) are real, A is an eigenvalue if and only if A is so. Therefore, the 
spectrum for matrices in 0(N,N) comes again as complex quadruples, or real pairs A, A -1 . 
These pairs are in general not degenerate though (an example for a matrix in 0(1, 1) with 
simple real eigenvalue pairs is the hyperbolic rotation). Similar as in fl3]), it is possible to 
define the groups 0(N,M) with maximal compact subgroup 0(N) x 0(M). Similar spectral 
properties as U(N, M) holds, but this will not be spelled out in detail. 



1.8 Class CII 

Class CII is also a subgroup of Q AIU , however, with a symmetry that requires tensorizing with 
another C 2 so that the matrices are of size AN x AN. The symmetry in U(2iV, 2N) is denoted 
by J cg> 1 instead of just J to highlight the difference. Then 

gan = | T G GL(4iV,C) \ T* J ®1T = J ®1, J ®PT J ®I = T) = SP(2iV,2iV) , 

and 



u cn = gen R TJ( 4i V) 



V 
w 



V,W E SP(2iV) \ = SP(2iV) x SP(2iV) . 



Obviously, Q CYl C U(2iV,2iV), but one also has 

D* g cu D = {T G GL(4iV, C)\T*I®IT = I®I, I* <g> IT I ® 1 = T } 

where D = 2^Q|) and this shows that D*Q 0ll D C U*(4A r ). It is not possible, however, to 
write Q cu as an intersection of U(2iV, 2iV) and U*(4A^) (which was possible in Class Dili). 

By the same reasoning as in Class Dili (see Section [L6l) the spectrum of T G SP(2iV, 2iV) 
comes in quadruples and real eigenvalues always have even geometric multiplicity. However, 
Tv = Xv here implies T J ® Iv = A J® I v. Class Dili has again a generalization SP(2A r , 2M) 
as in (J3J), but the obvious details are not given here. 



1.9 Class D 

Class D is again a subgroup of Q K of matrices of size N x N: 

g° = [T g GL(N, C) | T l T = 1 } = 0(N, C) . 

Note that this is also a *-group, namely T* G 0(N, C) if and only if T G 0(N, C). Clearly the 
maximal compact subgroup is U u = Q u D U(A r ) = 0(N). The relation T l T = 1 implies that 
the spectrum of T* is given by the ^-reflection A H- (A) _1 of the spectrum of T. Therefore the 
spectrum of a self-adjoint element T = T* G 0(N, C) is invariant under the map A i— > A -1 . 
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1.10 Class C 



Finally Class C is also a subgroup of Q A , albeit of size 2iV x 2N: 



{T G GL(2iV, C) | THT = I } = SP(2A, C) . 



This is also a *-group. The maximal compact subgroup is U° = Q° D U(2N) = SP(2A r ). The 
conclusions for the spectral theory of T £ SP(2A, C) are the same as in Class D. 

2 Lyapunov spectrum with random phase property 

Let be given a sequence (T n ) n >o of independently drawn and identically distributed random 
matrices in one of the above classical groups Q C GL(A r , C). For sake of simplicity, these 
distributions are supposed to be have fine moments (existence of third moment is sufficient for 
the second order perturbation theory later on). One is then interested in the growth in n of 
the norm of the products T n - ■ -Ti, or what is equivalent, the growth of their largest singular 
value. As one might expect, it turns out that this growth is exponential with a rate defining 
the top Lyapunov exponent: 



It is known |BL] that the limit converges almost surely to the same value and therefore an 
expectation can be introduced before the limit on the r.h.s.. Furthermore, instead of calculating 
an operator norm, an arbitrary vector v, say of unit length, may be inserted so that one actually 
only needs to calculate the length of vectors: 



Of course, also the growth rates 7 P of the other singular values are of interest. Together they 
constitute the Lyapunov spectrum. One convenient way to extract them is to use exterior prod- 
ucts A p Tn ■ ■ ■ T\ acting on the p-dimensional Grassmannian identified with the decomposable 
unit vectors in A p C iV where N is the size of the matrices in N |BL] . Then 



where V\,...,v p are unit vectors in which are at least orthogonal (further conditions will 
be imposed below when Q is not Q A ) and the norm is now calculated in A^C^. It is now very 
convenient to telescope the r.h.s. so that it becomes a Birkhoff sum. For that purpose let us 
introduce a sequence {v^ p \n)) n >Q of decomposable unit vectors in h?€, N by 



7i = 



lim - log(||T T1 ---T 1 ||) . 



7i 



lim -E log (||T n - -.T^H) . 



V lq = lim - E log (||A p T n ■■■T 1 v 1 A...Av p \\) 




APT n v(p\n - 1) 



W(0) = UiA...A«p. 



(4) 



v 



APT n v(p)(n - 1)|| 



v 
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Indeed, then one has 

p 1 n 

V 7(? = hm - E log(||A p T fc ^(fc-l)||) . 

* — ' n— s-oo ri ^ — * 
q=l k=l 

As the expectation over can immediately be taken in each summand (because v^(k — 1) 
is independent of J*,), one is therefore confronted with the calculation of Birkhoff sums of the 
random dynamical system (jl]) on the p-dimensional Grassmannian in which is a compact 
Riemannian manifold. It is known |BLj that this Markov process has a unique invariant and 
ergodic measure which we denote by v. In terms of expectations w.r.t. one then has the 
so-called Furstenberg formula: 

£7, = E r E, log(||A p T^)||) . (5) 

9=1 

The measure v is always Holder continuous and, if the distribution of T is absolutely continuous, 
then also v is absolutely continuous [BL] . Here we will be particularly interested in a situation 
where T n is of the form 

T n = R n e xp " , (6) 

where R n G Q is unitary, A is a real coupling constant and P n is in the Lie algebra of Q. If 
the distribution of the i? n and P n satisfies a certain coupling hypothesis formulated in |SS1] 
(of Hormander type), then the invariant measure v = v\ is in a weak sense and with errors 
of order A absolutely continuous w.r.t. to the invariant distribution on A^C^, and even equal 
to it for the so-called iV-orbital model of Wegner |SSlj . This invariant distribution is simply 
obtained by taking an element G A^C^ and randomly rotating it A P U with the unitary 
U drawn according to the Haar measure on the maximal compact subgroup U C Q. If Q = 
Q A = GL(iV, C) the measure on A^C^ is independent of the choice of the reference vector 
However, if Q is not equal to Q A = GL(JV, C) so that IA is smaller than IA A = XJ(N), then not the 
full Grassmannian is attained. Thus then also the choice of the reference vector becomes 
relevant. The guiding principle in determining the correct choice is that the vectors v q can be 
the eigenvectors of the a positive element in Q (which is (Tjv • • • Ti)*Tn ■ • • T\ in the application 
above). For example, if Q = U(A r , N), then this imposes that the v q are J-isotropic and pairwise 
J-orthogonal. Furthermore, if the spectral theory in Q leads to Kramers' degeneracy, then also 
the vectors v q have to respect that structure. The choice of the orthonormal basis will be 
discussed in detail in each case below. Apart from this (admittedly important) detail, we can 
now state the central hypothesis of this paper extending that of [RS2] . Strictly speaking it 
extends |RS2] to other symmetry classes, but is less general because no hyperbolic channels in 
the terminology of [RSI I IRS2j are considered here. 

Random phase property (RPP): The invariant measure v of the Markov process (jl]) gen- 
erated by a random sequence in Q is given by the invariant measure on the Grassmannian 
obtained by rotating an initial condition with unitaries drawn from the Haar measure on 
the maximal compact subgroup IA C Q. 
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We don't expect the RPP to hold exactly except in a model where it is imposed artificially. 
On the other hand, it may hold approximately in a perturbative situation described above and 
its implications can also be observed in numerical experiments on particular models, provided 
one restrict the attention to so-called elliptic channels (e.g. [RSI] ). As discussed in detail in 
[RS2], the RPP is strictly weaker than the maximal entropy Ansatz [Dorl IMPKj used for the 
derivation of the DMPK equations. A simple toy model for which the RPP holds exactly (and 
thus all the below as well) is obtained by requiring R n in (JHJ) to be Haar distributed onWc^ 
and then P to be distributed independently according to an arbitrary measure on the the Lie 
algebra of Q. Then it is obvious that the invariant measure v of the Markov process (|1J) is 
given by the Haar measure, namely the RPP holds. In general, provided the RPP holds, the 
calculation of the Lyapunov spectrum by Furstenberg's formula (J5j) becomes 

v 

Y,l q = EtE, logfllATA'Z/viA... Av p \\) . (7) 

q=l 

In a perturbative regime which is of interest for several applications, the evaluation of (J7J) is 
now an algebraic exercise. This was actually carried out in detail for the groups Q Am , Q CI and 
Q mu in |RS2j . Here we also consider the remaining cases which involves calculating averages 
not only over the unitary, but also the orthogonal and symplectic group. The basic tool will be 
the following perturbative formula in the case where T = Re xp with R unitary. Let v i, . . . , v p 
be ortho normal vectors in and set 7i p = X]q=i e q e q- Then 

log (\\A P Re xp A P U v\ A . . . Au p ||) 

A A 2 A 2 

= - Tr (U*Q Uir p ) + — Tr (U*S Uir p ) - — Tr (U*Q Un p U*QUn p ) + C(A 3 ) , (8) 

where 

Q = P + P* , S = P 2 + 2 P* P + (P*f . 
To prove this, let us begin by recalling the definition of the norm in A P( C N : 

l°g (\\A p Re xp A P U v%A... Av p \\) = \ log det p (v*(e XP U)*e XP Uv^ . . =1 ^ 
Next the identity log det p = Tr p log and an expansion in A shows equality to 

\ Tr p lo s (ip + A ( v i u *Q Uv j)i,j=i,..., P + y (<■:'-*' <■;),, , ; , + 0(\ 3 ) 

Now an expansion of the logarithm already proves (Ej)- Substituting (jHJ) into (JTj) one obtains 
a perturbative formula for the sum of the Lyapunov exponents. Taking the difference of two 
such sums then allows to deduce a perturbative formula for each Lyapunov exponent j p . This 
procedure and also the above proof of (jHJ) do not provide good error estimates on the dependence 
of the 0(X 3 ) term on the matrix size N. Actually, considerably improved bounds can be 
obtained using the Gram-Schmidt cocycle with values in the upper triangular matrices with 
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positive diagonal |RS2j . Further improvements in this direction will be provided elsewhere. 
Here the main focus is on the calculation of the contributions up to order 0(X 2 ) rather than 
the error estimate. Let us now come to the results in each of the 10 Cartan classes. Even 
though the final formulas are all quite compact, it ought to be added that the calculations 
leading to these formulas are quite tedious. We only provide the main intermediate steps. 



2.1 Lyapunov spectrum for Class A 

As this is the first case, let us treat it in detail. Due to Section 11.11 the group is GL(iV, C) 
with compact subgroup U(iV). As there are no specific symmetries, one may choose the unit 
vectors v q to be the standard basis vectors e q having only one non-vanishing entry in the qth 
component. Then replacing (jSJ) into ((7j) one sees that one needs to calculate second and fourth 
moments of the Haar measure on XJ(N). The formulas needed are listed in Lemma [TJ As 
Tr(7T p ) = p, one deduces with some care 

X>, A = E 



A pTr(Q) + A2 pTr(S) _ A2 (Np - p 2 )Tr(Q) 2 + (Np 2 - p)Tr(Q 2 



2N AN AN(N 2 



Now the identity Tr(S) = Tr(Q 2 ) allows to further simplify. Furthermore, taking the difference 
of Eq=i 1q and Eg=J 7?, one finds 



7 P A = A — E P Tr(Q) + A 4(jy2 _ ^ E P (Tr(Q 2 ) - - Tr(Q) 2 j + O(A') . (9) 

The main remarkable facts about this formula are the following. To order A, the Lyapunov 
spectrum is iV-fold degenerate. However, in second order perturbation theory this degeneracy is 
lifted and, furthermore, the spacings between the Lyapunov exponents are equal. This equidis- 
tance is observed in numerous numerical experiments and is here shown to be a consequence of 
the RPP. Finally, the Lyapunov exponent Tr V+i i is smallest in absolute value. Here [ . ] denotes 
the integer part. If iV is odd, then j fj^i , = if P is centered. However, this vanishing is lifted 

by higher order terms, just as a non-centered P would not lead to a non-vanishing Lyapunov 
exponent. Hence the vanishing is not symmetry enforced, as in Classes Dili and D with odd 
N. 



2.2 Lyapunov spectrum for Class AI 

Here all matrices are real as Q Al = GL(iV, K). Thus also the initial vectors v\, . . . , v p should be 
chosen real. Thus one can again choose v q = e q . Then one is lead to an analogous calculation 
as in Class A, but according to the RPP the averages now have to be taken over the orthogonal 
group instead of the unitary group. The corresponding formulas are collected in Lemma |2j 
After some algebra one finds 

jt 1 = A — EpTr(Q) + A 2 N + 1 ~ 2 P e (ti(Q 2 ) - — TrfQ) 2 ) + 0(\ 3 ) . 
The same comments as in Class A apply. 
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2.3 Lyapunov spectrum for Class All 

Here the matrices T = Re xp are in the group Q AU = U*(2N) and averages have to be taken over 
the symplectic group U AU = SP(2iV). However, some further care is needed with the choice of 
the Vi,...,v p . Indeed, the Lyapunov exponents are the scaling exponents of the singular values 
T n ■ ■ ■ Ti, thus the eigenvalues of the positive matrix (T n ■ ■ ■ Ti)*T n ■ ■ ■ T\ in the group U*(2iV). 
As explained in Section ll.3[ these eigenvalues are twice degenerate and the eigenvectors are 
pairs v, Iv. Thus also the Lyapunov spectrum is twice degenerate. If one chooses V\ = e\ and 
this leads to the largest singular value and thus largest Lyapunov exponent in (J7]), then one 
has to choose v 2 = Ie± so that indeed the second Lyapunov exponent is equal to the first one, 
hence establishing the degeneracy just explained. Arguing similarly for the other Lyapunov 
exponents, we therefore set 

v 2q -i = e q , v 2q = Iv 2q ~^l = I e q = e q+N . 

Then the projections 7r p = Y7 q =\ v q v q = i^pY are °f the f° rm 

j+5 



(o c) ' b = T,<w> c = E<«)*> ( 10 ) 

q=l q =l 



where p = 2j + 5 with 5 G {0, 1} and e' q denotes again the standard basis vectors of C^. Note 
that both b and c are N x N matrices. They satisfy be = cb = c as well as Tr(6) = j + 5 and 
Tr(c) = j. This implies 

r-K2j+8l 7T2j+S = 7T 2 j . 

Furthermore needed is the formula I*Q l IQ = I*QIQ = Q 2 . With these identities at hand, a 
careful use of Lemma [3] leads to 



v 



q =i 



A pTr(Q) ^ 2 pTr(S) 



AN 8N 

+ C(A 3 ) 



2 (2Np + 5- p 2 )Tr(Q) 2 + (2Np 2 - 2p - 2Np + 2NS)Ti(Q 2 
— A 



16iV(JV-l)(2JV+i; 



where 5 = |(1 — (— l) p ) as above. Now again Tr(S') = Tr(Q 2 ) allows to somewhat simplify. 
Then taking differences a further calculation shows 

,aii x 1 ^ ^ x2 2N + l-2p+(-l)r 



lt n = A E P Tr(Q) + A — r . / E P Ti(Q 2 ) Tr(Q) 2 + 0(\ 6 ) , 

Ip AN P w; 8(N -1)(2N + 1) P \ 2N w; J K 

where p = 1, . . . , 2N. Besides the comments made in Class A, let us note that this formula 
respects the double degeneracy of the Lyapunov spectrum resulting from Kramers' degeneracy. 



2.4 Lyapunov spectrum for Classes AIII 

The reflection symmetry of the spectrum of matrices in Q AIU = U(N, N) first of all implies 
that the Lyapunov spectrum satisfies 7p IH = — j^N+i-p- Hence it is sufficient to calculate the 
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q=l 

A 2 
16 



non-negative Lyapunov exponents. In principle, one can proceed as above, but it is important 
to choose the vectors v q J-isotropic because the Lyapunov spectrum studies the scaling of the 
singular values of products in U(iV, N), which are eigenvalues of positive matrices in U(iV, N) 
so that their eigenvectors are always isotropic. Therefore we choose 

v q = -j= (e q + e q+N ) . (11) 

It follows that 7T P = |(^) where d = Y7 q =i e 'q( e ' q )* is an N x N matrix built from the standard 
basis e' q in C^. Now the average over U AIU = U(iV) x U(A r ) can be calculated with Lemma d](i). 
Furthermore, the relation QJ + JQ = implies that Q = (°» g) is off-diagonal so that, in 
particular, Tr(Q) = 0. With the notations Il + = (Jjj) and II_ = (°°) one gets after carrying 
out the averages over the quadratic terms of U : 

= E P E y E iy [^^(Tr(Sn + )Tr( 7 r p n + ) +Tr(5n„)Tr(vr p n_)) 

(Tr(V* aWdW*a*Vd) + Tr (W* a* VdV*aWd))] + C(A 3 ) , 
so that Tr(7r p n±) = Tr(cZ) = p, Tr(5) = 2Tr(a*a) = Tr(Q 2 ) leads to 

q=l > ' 

Taking differences therefore shows 

N - v + - 

7 P Am = A 2 g ^ 2 EpTr(Q 2 ) + 0(X 3 ) , p=l,... t N, (12) 

which was already derived in |RS2] . One important feature of this formula is that the smallest 
non-negative Lyapunov exponent 7^ m is strictly positive. This does not hold if one considers 
random products in U(iV, M) with N < M. Indeed, then there are M — N vanishing Lyapunov 
exponents. With some care one can transpose the above calculation to show that the remainder 
of the Lyapunov spectrum is again given by (fT2~j) . but in the denominator iV 2 is replaced by 
NM. 

2.5 Lyapunov spectrum for Classes CI 

In Class CI £ CI ^ SP(2iV, E) and U C1 ^ U(iV). This was also dealt with in [RS2] and the result 



is 



,0, . ,2 N-p+1 



Then the other half of the Lyapunov spectrum is again given by 7°' = — 72jv +1 _ p - The proof is 
a mixture of Class AIII and Dili below. 
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2.6 Lyapunov spectrum for Class Dili 



In Class Dili the matrices are in the group G Am = XJ*(2N) n U(iV, N) = S0*(2iV) so that the 
Lyapunov spectrum has both the double Kramers' degeneracy and the symmetry around 0. For 
odd N, this implies that there must be a twice degenerate vanishing Lyapunov exponent. This 
vanishing is symmetry enforced and thus must also hold to all orders of perturbation theory. 
The implication of a vanishing Lyapunov exponent are discussed in Section |5j Now let us come 
to the perturbative calculation of the Lyapunov spectrum. For the choice of the orthonormal 
basis v q two criteria have to be satisfied. First of all, the eigenvectors have to be J-isotropic 



and pairwise J-orthogonal for different eigenvalues as in the case of Q f 
they should be chosen in pairs as in the case of Q An . Therefore, we set 



V 2 g-1 



1 

7! 



(e q + e q+N ) 



V2q 



I v 2q - 



1 

71 



and second of all, 



(13) 



Note that both these vectors are J-isotropic, but they are not J-orthogonal (which is not 
required because they correspond to eigenvectors of the same eigenvalue). However, each is 
J-orthogonal to all other basis vectors except for his pair partner. Then 7i p = J2q=i v q( v q)* * s 
precisely the same as in ffHJ]) . Now the average over the group U mu = U(N) in (J7J) is carried 
out using the formulas of Lemma HI Note that QJ + JQ = and I*QI = Q imply that 
Q = Q* = (°» ™) is off-diagonal with a 1 = —a, so that Lemma H](iii) applies. Also Tr(Q) = 0. 
Furthermore, n p is also of the form required in Lemma IU(ii) and (iii). With the notation II ± as 
in Lemma H] and writing again p = 2j + 5 with 5 G {0, 1}, 



E 

9=1 



A2 E 
-E P 



As one also has I* SI 
Thus 



s and j*n + 1 = n_ 



E 

9=1 



la 



A* 
16 



2p 
~N 



V 



it follows that Tr(n_<S') 
-2p + 5 



N(N 



E P Tr(g^ 



Tr(Q 2 

= Tr(n + 5) 
C(A 3 ) . 



0(A S 



^Tr(Q 2 



From this one deduces 



n -p + \ + \(-iy 

8N(N-1) 



E P Tr(g 2 



0(\ 3 



P 



N . 



Again the other half of the Lyapunov spectrum is given by 7°™ = — j^N+i-p- ^ s already 
pointed out above, the case of even iV was already treated in |RS2j . 

It is remarkable that the smallest non-negative Lyapunov exponents 7^r IH , 7^ and 7^™ for 
even N are strictly positive and satisfy 2^ IU = 7$ + C^iV" 1 ) and 2^ IU = 7^ in + C^A^" 1 ), 
the latter again in the case of even N. It reflects the change of the inverse localization length 
under breaking of time- reversal symmetry, see the discussion in |RS2] . 
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2.7 Lyapunov spectrum for Class BDI 

Because the spectrum of a self-adjoint matrix in Q Bm = 0(N,N) comes in pairs A, A" 1 , the 
Lyapunov spectrum always satisfies 7™* = —j^W-p- Hence it is again sufficient to calculate 
the first half of the spectrum. As 0(N, N) C U(iV, A), one has to choose the vectors v q as in 
(fTTj) . These vectors are already chosen to be real which would be the second requirement in 
the present case. Now both tt p and Q are as in Class AIII, tx p = |(^) where d = Y7 q =i e 'q( e 'q)* 
and Q = (°»q). The only supplementary property is that a is real. Also the calculation of the 
Lyapunov spectrum remains the same as in Class AIII, except that the average is now taken 
w.r.t. the Haar measure on 0(A) instead of U(A). As only the second moments enter and 
they are the same by Lemma and Lemma [H(i), also the final result of the calculation is the 
same: 

N - v + - 

7 p Bm = A 2 g ^ 2 EpTr(Q 2 ) + 0(A 3 ) , p=l,...,N. 
If one considers the case of 0(N, M), the same comments as at the end of Section I2T51 apply. 



2.8 Lyapunov spectrum for Class CII 

As already stressed in Section 11.81 elements of both groups Q cu and Q mu have Kramers' de- 
generacy and reflection symmetry, so that the Lyapunov spectrum is twice degenerate and 
symmetric around 0. However, the symmetry leading to Kramers' degeneracy is different and 
therefore instead of f|T3|) one sets for q — 1, . . . , N 



V2q-1 

From this one deduces 



-y= (e, + e q+2N ) 



V2q 



J ® I v 2q -i = -j= {e q+N - e q+ - iN ) . 



7r n 



e / 
/ e 



b 
c 



b 
-c 



with N x N matrices b and c defined in ( flOl) . Thus //// = c) • m particular, for p = 2j + 5 
one has 2Tr(e) = p and 2Tr(I fl f) = j. Furthermore, Q = (°*q) with Ial = a (note the sign 
which means that a is not a quaternion matrix, but xa is). Now from Lemma [5] one deduces 

v 



q=l 



A2 E 
-E P 



— (Tr(n + 5)Tr(n +7 r p ) + Tr(n„5)Tr(n_vr p )) 



4AT2 



Replacing thus shows 



E 

9=1 



7 = -r 



Tr(g 2 )(Tr(e) 2 + Tr(J7//)) + 0(\ 



E P Tr(g 2 ) + C(A 3 



A 2 


P 


- 1 


rp 2 




4 


4N 


4A 2 


V4 


2/ 



Therefore 



CII 



7n = 



A 2 2A-p + l 



-1? 



4A 2 



E P Tr(Q 2 ) + 0{X Z 



P 



2N . 
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2.9 Lyapunov spectrum for Class D 



The reflection symmetry of the spectrum of a positive matrix in Q v = 0(N, C) discussed in 
Section [L9l implies again that the Lyapunov spectrum always satisfies 7° = — 7^_ p . Therefore, 

if N is odd, there is one vanishing Lyapunov exponent 7^+1 = 0. This vanishing is hence 

2 

symmetry enforced. For the calculation of the Lyapunov spectrum, one starts again from (151) 
and (J7J). According the RPP the average now has to be calculated on the group U° = O(N), 
thus pends on the formulas in LemmaEJ Furthermore, one uses that Q = Q* G o(N, C) satisfies 
Q = —Q l = —Q so that, in particular, Tr(Q) = and Tr (Q t Q) = — Tr(Q 2 ). From this and 
Tr(S') = Tr(Q 2 ) one deduces after some algebra 



E 

q=l 



% = -7 



A_ 2 
4 



P 



p — p 



N N(N-l) 



X 2 



EpTr(Q^) - - 



Np — p 2 



4 N(N -1)(N + 2) 



EpTr(g) 2 + 0(\ s 



which implies 



7 P = -r 



A 2 N + 1 - 2p 
1 N(N - 1) 



E, 



Tr(Q 2 



1 



N + 2 



+ 0(X 3 ) 



P 



1 N . 



Hence in Class D there is never a term of order 0(A), and the perturbative formula indeed 
respects 7^+1 = 0. 



2.10 Lyapunov spectrum for Class C 



l2N-pi as is even, this 



Similar as in Class D the Lyapunov spectrum satisfies ^ = 
symmetry enforces no vanishing of a Lyapunov exponent. As preparation for the calculation of 
the Lyapunov spectrum, let us note that Q = Q* G sp(2A^, C) satisfies Q = —l*Q l l = —I*QI 
and thus Tr(Q) = and Tr (FQflQ) = -Tr(Q 2 ). By the RPP, the average in © is calculated 
over the group IA° = SP(2A^), based Lemma [31 One finds, first for p < N, 



E 

q=l 



7 C 



A 2 

— E P 



2^ 



P 



2N{2N 



Tr(Q 2 



+ C(A 3 



and therefore 



7 D = -r 



A 2 N+l-p 



E P Tr(Q 5 



C(A 3 



P 



N 



4 JV(2JV + 1) 

The second (negative) half of the Lyapunov spectrum is given by the reflection stated above. 



3 Symmetries of one-particle Hamiltonians 

This section briefly reviews the ten-fold way of Altland and Zirnbauer [AZj . albeit with a 
classification of the symmetries in the spirit of |SRFL| IRSFLj . In this classification the one- 
particle (first quantized) Hamiltonian H acting on a complex Hilbert space can have three 
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even TRS) 


J*HJ = 


H, 


J 2 = 


1, 


(odd TRS) 


PHI = 


H, 


I 2 = 


-1 


even PHS) 


K*HK -- 


-H, 


K 2 = 


1, 


(odd PHS) 


L*HL = 


-H, 


L 2 = 


-1 


(SLS) 


M* H M 


= -H, 


M 2 = 


±1 



basic symmetries, the time reversal symmetry (TRS), the particle-hole symmetry (PHS) and 
the sublattice or chiral symmetry (SLS). The first two of these symmetries can be either even 
or odd. These symmetries are implemented by real unitaries J, I, K, L and M in the following 
manner: 

(14) 
(15) 
(16) 
(17) 
(18) 

The notational choice of letters J, /, etc., is not yet connected to the symmetries as defined in 
(JTJ), but will in some cases be so in our examples later on. For instance, it is possible to have 
J = 1 so that the even PHS just reads H = H. 

Let us add a few comments. A given Hamiltonian H can have only one TRS and one PHS 
symmetry. If one TRS and one PHS is present, then one also has a SLS symmetry. For example, 
if (fTol) and (TIB"]) hold, then also (1 18)) with M = K*I. This also imposes the sign of the SLS which 
is why one does not keep track of it. The only case in which this sign may seem interesting is 
the case where there is only SLS. But in this situation the reality of M is irrelevant and then 
M may be changed to i M which alters the sign. Counting all possibilities to combine the three 
symmetries, one obtains 10 classes which are listed in Table 1 (1 with no symmetry, 5 with 
one symmetry and 4 with three symmetries). The table also contains the Cartan classification 
associated to each such H Class by Altland and Zirnbauer |AZj in the following manner. Each 
Cartan class is specified by an involution on a compact Lie algebra. The —1 eigenspace of 
the involution multiplied with i then produces a set of self-adjoint matrices which forms the 
corresponding H Class of Hamiltonians. See also |HHZ| ISRFLl IRSFLl ISCRl lAKj for further 
insights on this connection. 



4 Quasi-one-dimensional model operators 

In this section we construct for every H Class in Table 1 a quasi-one- dimensional model operator 
of Dirac type. If the symmetries are organized using the invariance under representations of 
the real and complex Clifford algebras [HHZ[ IKitl ISCRl lAKj , then there is a nice constructive 
way to obtain these model operators. Even though this might be quite illuminating, it is not 
the focus of the present paper where we only restrict ourselves at providing one example for 
each given class, similar as in [BMSA, iBFCMl ITBFMj . Our model operator of Dirac type will 
have a random perturbation and the study of its fundamental solutions (transfer matrices) then 
leads to one of the classes of random products described in Section [2J The model operators 
act on the Hilbert space L 2 (IR, C 2N ) or, if necessary, on L 2 (IR, C 2N ) ® C 2 , and are of one of the 
three following forms: 

Hj = Jid - AV, Hi = Id - XV, H K = Kid - XV , (19) 
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H Class 


TRS 


PHS 


SLS 


Q Class 


g 


U = 


H 


A 











AIII 


U(N, M) 


U(JV)xU(M) 


1 


AI 


+1 








CI 


SP(2N, E) 


U(iV) 


1 


All 


-1 








Dili 


SO*(2A0 


U(N) 


2 


AIII 


o 


o 


1 


A 


GL(iV CI 


U(N) 


1(21 


BDI 


+ 1 


+1 


1 


AI 


GL(N, E) 


0(N) 


1(2) 


CII 


-1 


-1 


1 


All 


U*(2N) 


SP(2N) 


2(4) 


D 





+1 





BDI 


0(N, M) 


0(JV)xO(M) 


1 


C 





-1 





CII 


SP(2N,2M) 


SP(2iV)xSP(2M) 


1 


Dili 


-1 


+1 


1 


D 


0(iV,C) ' 


0(N) 


2 


CI 


+1 


-1 


1 


C 


SP(2A^, C) 


SP(2N) 


2 



Table 1: Table with symmetries taken from jTBFM} ISRFL] together with the transfer matrix 
group Q for quasi-one-dimensional operators H and its maximal compact subgroup IA. Further- 
more, fi is the symmetry imposed multiplicity of the Lyapunov spectrum given by the spectral 
multiplicity of self-adjoint elements of Q. The numbers in parenthesis give the multiplicity to 
lowest order in perturbation theory for centered perturbations. 



Here J, / and K are the matrices given in ([T]) acting on C 2N and the supplementary C 2 allows 
to implement a further symmetry (this is only necessary for H Classes C, CI and CII because 
an odd PHS cannot be implemented in any of the three model Hamiltonians) . Furthermore, d 
is the derivative on L 2 (M, C 2N ) and V = V* is a matrix- valued potential. It could be a matrix- 
valued function, but for sake of concreteness let us rather only consider a potential given by a 
sum of Dirac peaks: 

V(x) = J2 V nS(x-n) . (20) 

The V n are independent and identically distributed random matrices. The formal definition 
of the singular potential is via boundary conditions, namely either ip(n+) = e XlJVn ip(n—) 
or %p(n+) = e XIVn ip(n—) or ip(n+) = e XlKVn ip(n—), depending on which model in (TH?]) is 
considered. These formulas combined with Weyl-Titchmarch theory allow to show that Hj, Hj 
and Hk are well-defined self-adjoint operators (see |SS2j for details). 

In H classes A,C and D it is possible to consider Hq = Gid — A V acting on L 2 (R, C N+M ) 
where G is defined in ([3]) and < M. This models a system with a different number of left and 
right movers and gives an effective description of chiral edge states of quantum Hall systems 
(conventional, spin and thermal quantum Hall effect respectively), which is reflected by the fact 
that the Lyapunov spectrum contains at least M — N vanishing Lyapunov exponents. These 
systems will be discussed briefly in Section 

For each real energy E, the fundamental solutions T E (x, y) of the Schrodinger equation 
HT E (x,y) = ET E (x,y) which are right continuous in both x and y and satisfy the initial 
condition T E (y,y) = 1 are also called the transfer matrices of the system. By the above, they 
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are given by 

T E (n-,l-) = T n ---T 1 , (21) 
where the matrices T n for three cases in (TT9l) are respectively given by 

rp _ lEJ XlJVn rp _ E I XIV„ T _ tE K XlKVn 

± n CO , ± n CO , J. n CO 

These matrices are of the form if one sets R n = e lEJ and P n = i JV n in the first case, 
and correspondingly in the other cases. Now by construction, the transfer matrices are in the 
following three groups: 

Tj = U{N,N), Ti = C*XJ(N,N)C , Tk = C\J(N, N) C* , (22) 

where C is the Cayley transform defined in ()2]). For H Class C, CI and CI I one has to replace 
U(N,N) by U(2iV, 2N) due to the supplementary C 2 . Thus it looks like the transfer matrices 
are always in the Cartan Class AIII, but actually the game in the following subsections is 
to show how the symmetries of the Hamiltonian (TRS, PHS, SLS) lead to subgroups of Tj, 
Ti or Tk which then actually cover all Cartan classes. Table 1 collects the results of this 
correspondence. 

4.1 Standard unitary class (H Class A) 

Let us choose the model Hj. As there are no further symmetries, the transfer matrix group 
Tj is equal to Q Aul = U(N, N) and therefore the results of Section 12.41 can be directly applied 
to study the associated Lyapunov spectrum. 

4.2 Standard orthogonal class (H Class AI) 

Let us again choose the model Hj and then implement the even TRS by K*HjK = Hj. Then 

K*%K = K*e~ lEJ e~ xuW K = e lEJ e XlJK *^ K = T n , 

shows that Tf 1 = Q 01 = C SP(2A r , R) C* and Section [231 applies. Alternatively, one can choose 
H! and the implement the even TRS by = Then 7} AI = SP(2iV, R). 



4.3 Standard symplectic class (H Class All) 

Starting again from the model Hj, the odd TRS is now PHjI = Hj. By a similar calculation 
as above one find that T/ n = G mu = CSO*(2N) C*. Thus Section ES applies. 



4.4 Chiral unitary class (H Class AIII) 



When there is a chiral symmetry, it is good to work with the model operator Hj and implement 
the SLS by J*HjJ = —Hj. The reason is that this implies J*T n J = T n , which combined with 
/ shows that T n is in the group 



T*IT, 



Tf 



A 






(A- 1 )* 



A e GL(N, C) 



(23) 
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Therefore 7^ Am ^ g A and the results of Section |2~T1 for Class A apply. The formula in Section [27T1 
gives the Lyapunov spectrum of the upper component in 7} Ain . The exponents of the lower 
component are then given exactly by the negative of the upper component. In this manner 
one recovers the symmetry of the Lyapunov spectrum which has to hold for the group 7} Am C 
XJ(N,N). It is worth mentioning the following implication of the comments at the end of 
Section 12.11 If the potentials V n are centered with an distribution that is also even, then 
already the Lyapunov spectrum of the upper component given in (Q is symmetric around up 
to terms of order 0(X A ). Therefore, the Lyapunov spectrum in 7} Am has a double degeneracy up 
to terms of order 0(X 4 ). If, moreover, N is odd, then there are two Lyapunov exponents that 
vanish to second order perturbation theory and the first non-vanishing contribution appears 
only at order 0(X 4 ). This leads to a very large localization length in such systems. The same 
comments apply to the other two chiral classes. 



4.5 Chiral orthogonal class (H Class BDI) 



Everything of the last section transposes, but, moreover, one has the even TRS Hi = Hj. 
Combined one also obtains an even PHS J* Hi J = —Hj. The transfer matrices are now, on top 
of being in 7} Am , real as in Section 14.21 and therefore 



r/ DI = 

Thus the results of Section 



A 




A e GL(N,\ 





apply, as well as the comments made in Section 14.41 



4.6 Chiral symplectic class (H Class CII) 

As in the other two chiral classes, the Hamiltonian is Hj with SLS J*HjJ = —Hi, but now 
it acts on L 2 (M, C 2N ) (g) C 2 and the supplementary C 2 is used to implement the odd TRS by 
1 <g) I* Hi 1 <gi I = Hi (note that here the / £g> 1 in the definition of Hi and 1 <g> / in the 
odd TRS are in different gradings). As usual the two symmetries together imply another one, 
namely an odd PHS J ® I* Hi J ® I = —Hi. One deduces that the transfer matrices satisfy 
1 <g> I* % 1 <g> I = T n so that 



A e u*(2iv) } g 



ran / (A 

"Ho (A" 1 )' 

Thus the results of Section 12.31 apply, which lead in particular to Kramer's degeneracy, which 
leads even to a fourfold degeneracy up to order C(A 4 ) if the distribution of V n is even. 



4.7 Even PHS without TRS (H Class D) 

Here let us choose the model Hj and implement the even PHS simply by Hj = —Hj. This 
implies for the transfer matrices that T n = T n , which leads to the following subgroup of Tf'- 

Tf = {T e GL(2iV, C) | T* JT = J , T = T } = 0(N, N) = £ BDI . 

Thus the formulas of Section 12.71 apply. 
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4.8 Odd PHS class without TRS (H Class C) 

The odd PHS cannot be implemented in any of the three model Hamiltonians, unless they are 
tensorized with an extra C 2 . Having done that, let us choose the model Hj. Then the odd 
PHS is given by J <g> I*TTj J ® I = -Hj. From this one deduces J <g> I* % J ® I = T n so that 
the transfer matrix group is 

Tf = {T G GL(4N,C) | T*J ® IT = J®1 , J ® I*T J ® I = T } = Q cu . 

Thus the Lyapunov spectrum can be calculated as in Section 12.81 

4.9 Even PHS with odd TRS (H Class Dili) 

In Class Dili let us use the model operator Hk and implement the odd TRS by I*HkI = Hk 
and the even PHS by K*HkK = —Hk- Then one readily checks that the transfer matrices are 
in the following group: 

7£ m = {T e GL(2iV,C) | T*KT = K , I*T I = T , K*T K — T} . 

Writing out the latter relations implies (in particular, first deduce and then use J*TJ = T) 

A e 0(N,C) | g D . 

Thus the results of Section 12.91 apply. In particular, the Lyapunov spectrum of Hk is always 
twice degenerate (because the spectrum in 0(N, C) is symmetric) and there is always a twice 
degenerate vanishing Lyapunov exponent if iV is odd. 

4.10 Odd PHS and even TRS (H Class CI) 

Again the odd PHS enforces the use of a supplementary grading, so the Hilbert space is 
L 2 (IR, C 2N ) <S> C 2 . Let us choose the model operator Hk and implement the even TRS by 
I (g) IHk~I (g) I = H K and the odd PHS K ® I*Hk~K ® I = —H K . Then also the SLS 
J ® 1 Hk J <8> 1 = — Hk holds. Next one verifies that the transfer matrices are in the following 
group: 

T^ 1 = {T e GL(4iV,C) \T*K®1T = K®1, K ® I*T K ® I = T , J ®1T J ®1 = T) . 
After some algebra, the latter relations show 

A G SP(2A^,C) | = Q c . 

Hence the results of Section 12.101 apply. In particular, the Lyapunov spectrum of Hk is always 
twice degenerate (because that in SP(2A r , C) is symmetric). 



«tdih / (A \ 
' K ~ \ \ (A- 1 )* J 



r c: _ r A \ 
' K - 11 (A' 1 )* J 
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5 Symmetry enforced delocalization 



The paper |SS2j studies the model Hamiltonian Hi of ( Tl9l) in the symmetry class All and, in 
particular, the case where the channel number N is odd. As explained in Section 14.31 there are 
at least two vanishing Lyapunov exponents in this case, and if a certain coupling hypothesis 
holds, there are exactly two vanishing Lyapunov exponents, just as in the RPP. By Kotani 
theory the two vanishing Lyapunov exponents imply that the spectrum of H has almost surely 
an absolutely continuous component of multiplicity 2. Furthermore, it was possible to show in 
[SS2J that almost surely there is no singular spectrum provided the distribution of the V n is 
absolutely continuous. All these results actually transpose to H Class Dili in the case where 
N is odd. No details of proof are provided here because there is really no essential difference. 

Finally let us discuss systems with a different number of left and right movers which are 
relevant for the modeling of the edge modes in various quantum Hall systems. 

Hq = Gid + V , 

acting on L 2 (1R, C N+M ) where N < M are integers. If H satisfies a TRS or a SLS, then N = M 
because otherwise the kinetic part is not isospectral to its negative. But in the remaining H 
Classes A,C and D it is permitted to have M > N, see Table 1. It follows that the transfer 
matrices are in the group U(iV, M). For the H classes C and D they are actually in a subgroups 
0(N,M) and SP(N/2, M/2), where N and M have to be even in the latter case. Now the 
spectral theory of a self-adjoints T e U(iV, M) shows that T has at least M — N eigenvalues 1, 
see Section ll.4[ Thus Hq must have M — N vanishing Lyapunov exponents. Under a coupling 
hypothesis on the random potential there are no other vanishing Lyapunov exponents. Again 
it follows from Kotani theory alluded to above that these Hamiltonians have almost surely an 
absolutely continuous spectrum of multiplicity M — N . Again one can show that this absolutely 
continuous spectrum is almost surely pure if the distribution of the V n is absolutely continuous. 

Acknowledgements: We are thankful for financial support by the DFG. The basic idea to 
write this paper resulted from discussions at an SFB/TR 12 meeting in Gdansk in 2009. We 
all thank M. Zirnbauer for the invitation. H. S.-B. also thanks the Mittag-Leffler Institute for 
a quiet office space which finally allowed to finish this work. 

A Moments for the Haar measures over compact groups 

Section IA.1I of this appendix collects formulas for certain averages over the compact classical 
groups that are used in the main part of this paper. These formulas can be derived rather easily 
from results scattered over the literature, and those involving averages over the unitary group 
were already listed in |RS2j . Nevertheless, it is emphasized in the remaining sections of the 
appendix that all these formulas follow from a useful calculus involving the so-called Weingarten 
function [Weil ICSl ICSt] which may be viewed as a systematic approach to invariant theoretic 
arguments in the physics literature on averages, as in Section VIII of |AZj . To illustrate how 
invariant theory is applied, Section IA.2I deals explicitly with the fourth moments over the 
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symplectic group (averages of a product of four matrix entries) and provides the corresponding 
value of the Weingarten function. This allows in Section IA.3I to deduce Lemma |3]^ii) which is 
actually the most complicated of all formulas in Section IA.1I All others follow similarly from 
prior results on the orthogonal group [CSj . Higher moments can also be calculated explicitly, 
provided that the matrix size is sufficiently large compared to the order of the moment (that is, 
the number of matrix entries that are being considered). Section [A. 4l presents the general version 
of the Weingarten integration formula for the unitary, orthogonal, and symplectic groups. These 
results are not really used in this paper, but they are definitely of independent interest and are 
provided for further reference. 

A.l Collection of formulas used in Section [2] 

Lemma 1 Let A, B,C,D G Mai(N x N,C). The following holds for averages over U(N): 

(i) (Tr(U*AUB)) 

(ii) ( Tr(UAUB) } 

(iii) ( Tr(U*AUBU*CUD) ) 



(iv) ( Ti^AUBU'CUD) ) 



(v) (Tr (WAUBU'CUD)) 



All formulas remain valid if all U's are replaced by their complex conjugates, e.g. 

( Tr(U t AUBU t CUD) ) = (Tr(U*AUBU*CUD)) . 

Lemma 2 Let A, B,C,D G Mat(A r x N,C). The following holds for averages over O(N): 
(i) 

(Tr (O* AO B)) = — Tr(A)Tr(5) , (Tr(OAOB)) = — Tr(AB*) . 



lTr(A)Tr(5), 

^ Tr(AB*) , (Tr^AUB)) = 0, 

r Tr(A)Tr(C)Tr (BD) + Tv(AC)Tr(B)Tr(D) 
Tr(AC)Tr(BD) + Tr(A)Tr(5)Tr(C)Tr(D) 



N 2 - 1 
1 



AT(N 2 - 1 

N ^ - Tt(A)Tt{C)Tt(BD) + Tr(AC t )Tr(BD 



N(N 2 - 1) 
1 



AT2 



N(N 2 - 1) 



Tr(AC )Tr(BD) + Tt(A)Tt(C)Tt(BD ) 
Ti{AC t )Ti(BD t ) + Tr (AC )Tr(B)Tr(D) 

Tr(AC)Tr(BD t ) + Tr(AC t )Ti(B)Tr(D) 
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(ii) 



TriO'AOBO'COD)} 
N + l 
N(N -1)(N + 2) 
1 



Ti(A)Tr(C)Tr(BD) + Tr(AC")Tr(££>') + Tr(AC)Tr(B)Tr(D) 
Tr(A)Tr(C)Tr(BD t ) + Ti(A)Ti(C)Ti(B)Ti(D) + Tr(AC")Tr(.B.D) 



N(N -1)(N + 2) 

+ Ti(AC t )Ti(B)Ti(D) + Tr(AC)Tr(BD) + Tr(AC)Tr(BD 



(iii) 



( TiiO'AO'BOCOD) ) 
N + l 
N(N -1)(N + 2) 
1 



Tr(A*5C rf D) + Tt(DCBA) + Tr(A(7)Tr( J B)Tr(D) 
Ti:(A*BC*D*) + Tr(A t BC t )Tr(D) + Tr(OBAD* 



iV(iV- l)(JV + 2) 

+ Tr(CBA)Tr(D) + Tr(CAD t )Tr(B) + Tr(CAD t )Tr(B) 



Lemma 3 Let A, B,C,D G Mat(2iV x 2N, C). The following holds for averages over SP(2iV)-' 
(i) One has 

(Tr(U*AUB)} = ^Tt(A)Tt(B), (Tt^AUB)) = ± Tr(IA) Tr(J* B) , (24) 
and 

(Tt(UAUB)} = (Tt(U*AU*B)) = ■^—Tr(AIB t I) , (Tr (UAUB)) = — Tr(AB*) . 



(ii) 



Tr(U*AUBU*CUD)) 
1 



L (2iV - 1)Tt(A)Tt(C)Tt(BD) + Tt(A)Tt(C)Tt(I* B I D) 

-Tr(A)Tr((7)Tr(£)Tr(L>) - Tr(/M*/C)Tr(fi£)) + Tr(/*v4*/C)Tr(£?)Tr(.D) 
-Tr(AC)Tr(RD) - (2JV - l)Tr (PA* I C)Tr(I* B l ID) 

—Tv(AC)Tr(I*B t ID) + (2N - l)Tr(A(7)Tr(5)Tr(L>) . 

Let us close this section by stating two results that follow from the above. Recall from 
Section [1.51 the special form of U C1 = U mi1 C U(2iV) which is is isomorphic to XJ(N). The 
following lemma presents no general formula for second and fourth moments of this group, but 
only those moments which are needed in Section 12.61 The proof of the lemma is a calculation 
based on the identities in Lemma [TJ 
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Lemma 4 Let A,B e Mat(2iV x 2N,C). Then the averages over W l = U mu lead to the 
following: 

(i) With n + = (J J) and n_ = (J J) 

( Tr(U* AUB) ) = — Tr (II+A) Tr (n+B) + — Tr (n_A) Tr (]!_£?) . 

(ii) Let A and B be of the form 

. ( a\ _ AO 

where all entries are of size N x N. Then 

(Tr(U*AUBU*AUB) ) = Z Tr ( bc< )) + ^g^i Z ^(b)Tr(c)) 

{ 1 J ' iV(iV 2 - 1) 

(iii) //j moreover, a 1 = —a, then 

Tr(A 2 )(Tr(6)Tr(c) - Tr(6c')) 



Tt(U*AUBU*AUB)} 



N(N-l) 



Recall from Section [1751 the special form of U cu = SP(2iV) x SP(2A r ). The following lemma 
presents no general formula for second and fourth moments of this group, but only those 
moments which are needed in Section 12.81 The proof of the lemma is a somewhat tedious 
calculation based on the identities in Lemma Etl). 

Lemma 5 Let A, B e Mat(4A r x 4iV, C). Then the averages overU 01 lead to the following: 

(i) With n + = ( Q °) and Il_ = ( °) as above, 

(Tt(U*AUB)) = Tr (U + A) Tr (Tl + B) + ± Tr (n_A) Tr (n_S) . 

(ii) Lei A and £? 6e o/ i/ie /orm 



where all entries are of size 2N x 2N and furthermore Ial = a and B* = B f = B. Then 

TrM 2 )(Tr(e) 2 + Tr(r//n) 
(Tr(U*AUBU*AUB)) = — — n v } — v LJJl m 
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A. 2 The symplectic Weingarten function for fourth moments 

As already pointed out, the formulas of Section IA.1I can be deduced from |Colt ICS] except 
for averages involving the symplectic group. In the symplectic case a sign factor is missing in 
|CS] . This was corrected in |CStj . but no explicit value of the Weingarten function was given. 
This section treats the fourth moment needed for the proof of Lemma E](ii) in detail. This also 
indicates how in principle the general results stated in Section IA.4I can be proved. 

Let ej, j = 1, . . . , 2N, denote the standard basis of the complex vector space V = C 2N . For 
any linear map T on V, the matrix entries are denoted by = e*Tej. Associated to / defined 
in (JTJ), let us introduce on V the skew symmetric bilinear form a(v,w) = w'/u; where v l = v* 
denotes the transpose of v £ V. Then for i, j = 1, . . . , N, 

Ti+N,j = — a (ej,Tej), T i+N j +N = — a(ei,Tej + N) , 
Tij = a(ei + N,Tej) , T^+n = a(ei + N,Tej + N) ■ 

Thus, for . . . ,z' 4 , j 4 £ {1,...,N} and . . . , a 4 ,/3 4 £ {0, 1}, 

/ dU t/j 1+Q!1 iVj 1 + / 3iJV Ui 2+ce2 N,j 2 +P2N Ui 3 +a 3 N,j 3 +l3 3 N £4 1 +«4A r j 4 +/34A r (^) 

JSP(2N) 

= (_ 1 p+a 2 +«3+«4 a ( eh+{1 _ ai)N o . . . (gj e i4+( i_ a4)7V , / dU Ue jl+hN ® . . . ® t/e i4+/347V J , 

V JSP(2V) / 

where a is extended to V® 4 via 

a(v\ ®v 2 ®vz® v 4 , wi <g> m; 2 g) u; 3 g) 104) = a(t»i, ioi) a(w 2 , w 2 ) a(w 3 , u; 3 ) a(w 4 , w 4 ) . 



Note that this yields a symmetric bilinear form on V® . By the definition of Haar measure, the 
integral 

/ dUU e jl+/ 3 lN ® . . . (8) U e h+ p iN (26) 

</SP(27V) 

is a SP(2A r )-invariant in V® 4 , specifically, it is fixed under the action of U £ SP(2A^) on V® 4 
which on decomposable tensors is given via U ~(vi®v 2 ®v 3 ®v i ) = {Uvi®Uv2®Uv 3 ®Uv i ). If N > 
2, then by the symplectic case of Weyl's First Fundamental Theorem for tensor invariants |GW| 
Thm 5.3.3] and the results of Section 3.4 in [StoJ, a basis of the (three-dimensional) subspace 
of invariant tensors in V® 4 can be described as follows. Write pi = {{1, 2}, {3, 4}}, p 2 = 
{{1, 3}, {2, 4}} and p 3 = {{1,4}, {2, 3}}. These are the pair partitions of the set {1, 2, 3, 4}. To 
control the sign factors that arise from the skew symmetry of the form a on V, let us keep track of 
the natural order on the individual blocks by using the following notation. For m £ {pi, p 2 , p 3 }, 
let us write m = { (mi (m) , ri\ (m) ) , (m 2 (m) , n 2 (m) ) } , where {mi (m) , ri\ (m) } , {m 2 (m) , n 2 (m) } are 
the blocks of m, and m^m) < ni(m), m 2 (m) < n 2 (m). To pi is now associated the tensor 

N 1 

ft = Opi = z2 z2 e r?i+ ei 7v ® e r , 1+ (i_ ei)A f(-l) ei ® e m+e2N ® e m+{ i- e2)N {-iy 2 . 

r/l,»?2=l £1,62=0 
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Here, the maps r i— > 7] r and r h-> e r . are constant on the block {mi(pi), rii(pi)} as well as on 
{m 2 (pi), n 2 (pi)}, and to the m component of a block corresponds a vector of the form e nr+er N, 
while to the n component corresponds a vector of the form (— l) er e r?r+ (i_ er )Ar. Analogously, 

N 1 

#2 = 9p 2 = ^2 Si+eiJV ® e^ 2+e2 jv ® e m + (i_ £l )7v(-l) £1 <g> e^ + (i„ e2 ) A r(-l) e2 , 

??1,7?2=1 £1,62=0 

AT 1 

#3 = #p 3 = S e r?i+,iiv®e^ +e27 v®e^ +( i_ £2 ) 7 v(-l) e2 ®e J?1+{1 _ ei) jv(-l) ei . 

?71,?72=1 £1,62=0 

Then # 2 , #3} is a basis of the space of symplectic invariants in and it is a straightforward 
(tedious) computation to verify that the Gram matrix of the form a with respect to this basis 
is given by 

AN 2 2N -2N 
2N AN 2 2N 
-2N 2N AN 2 

As long as N > 2, this matrix is invertible with inverse 



1 



2N-1 -1 1 

mN-wN+vy- 1 2N 1 ~ 1 ~^ r _ 1 I ■ < 27 > 

By definition, the (z,j)-entry of this inverse is the value Wg SP ( 2A r)(pj, Pj) of the symplectic 
Weingarten function. In Section \A.A\ the general Weingarten function Wg SP ( 27V ) is defined by 
the matrix entries of the inverse of the Gram matrix of a w.r.t. a basis of the invariants, which 
in turn is in bijection with the pair partitions (of 2k points if 2/cth moments are considered). 
In order to evaluate ( I25p . one now expresses the invariant ( 12 6 p in terms of the basis {#1, # 2 , O3}. 
Noting that a(e il+aiN ®. . .®e iA+aAN , 6 p .) vanishes unless i mr ( Pj ) = V( Pj ) and a mr{Pj) = l-a nr{Pj) 
for r = 1,2, one obtains 

/ dU U^+^NJ^^N Ui 2+a2 NJ 2 +/32N Ui 3+a3 N,j 3 +/3 3 N Ui 4+a4 N,j 4 +l34N (28) 

JSP(2N) 

= ^2 Wg SP ( 2A r) (m, n) (— lj^iW+^w+^iW+^sW 

m,ne{pi,p2,p3} 

^('mi(m)i V(ra)) ^(^mi(ni)) 1 ^ni(m)) $v"m,2{xti)i ^ri2(m)) ^m2(m)) 1 ^n2(m)) 

) Jn 2 (n) ) 5(/3 m2 ( n ), 1 - /3 n2 (n)) • 
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A. 3 Proof of Lemma Mil) 

Let us begin by writing out the l.h.s. explicitly: 
(Tr(I£/*i AUBltflCUD)) 

N 1 

= ^ ] ^ ^ Iii+aiN,ii+(l-ai)N Ii2 + {l-a2)N,i2+a2N ^i2+a2N,i3+asN Bi 4 ^. a4 j\f t i s ^. a5 ]y 

il,«2,. ■■1*8 = 1 ai,ce2,...,ag=0 

Iis+a 5 N,i5+(l-a 5 )N I i(i + (l-a<i)N ,is+aaN Ci e> + a6 Nj 7 + a7 N -Dj g +a g Af,ji+ai N 
(Ui 2 +(l-oi2)N,ii+(l-a 1 )N ^3+03^.14+04^ ^6+(l-a6)A r ,i5 + (l-05)A r ^7+07^8+08^) • 

The next step is to expand this average ( . ) = J SP ^ 2N ^dU according to (I28I) . Then one 
interchanges the sum over pairs of pair partitions with the sums over the i and a indices and 
then determines the contribution of each pair (pi,pj). Let us illustrate this procedure for the 
pair (pi, p 2 ). One reads off from (128]) that for a choice of i and a indices to give a non- vanishing 
contribution, one must have that i 3 = ii, a 3 = a 2 , %-j = a 7 = a 6 , i 5 = ii, a 5 = 1 — a 1; i 8 = 
i 4 , as = 1 — «4. Furthermore, the sign factor corresponding to each choice of i and a indices is 

f]\Q-— <*2)+(l— a6)+(l— ai)+«4 _ /_-^\ai+Q2+Q!4+(l— cte) 

Consequently, the coefficient of Wg SP ( 2 j\n(Pi) P2) in the sum over pairs of pair partitions is 

E, 1 \ai+a2+a4+(l~ae) j j A r 

l — 1 ) - , ii+oiA r ,ii + (l-oi)Af 1 i2+(l-a 2 )N ,i 2 +a 2 N n -i2+a 2 N,i2+a 2 N n i 4 +a 4 N,i 1 +(l-a 1 )N 

i.o: 

Iil+(l-ai)N,h+aiN Ii(;+(l-a(;)N,i(j+a(;N C'i6+Q!e-'V,i6+a!6-'V -Di4+(l-a4)7V,ii+ai N • 

Observing that Ii+ a N,i+(i-a)N = (— l) 1 "" and h+(i- a )N,i+aN — (~ 1)°; one can group the factors 
as follows: 

7 4i2+ct2-?V,*2+a2 7V , 
Ci 6 +a 6 N,i 6 +a 6 N j 

-Si 4 +a4A r ,U+(l-ai)A r ^ii+(l-ai)JV,ii+aiiV -Di4+(l-a4)A r ,n+aiA r ( — l)"" 



= -Si 4 +a4A r ,U+(l-ai)A r -^ii + (l-ai)A r ,u+aiA r -D't4 + (l-a4)A r ,«i+aiA r ^4+(l-04)A r ,«4+«4A r ) 

Al+Q!liV,il + (l-Ql)iV ( _ I)" 1 = _ 1 ? 

^i2 + (l-a2)A r ,«2+a2A r ( — l)" 2 = 1 5 

Ii 6 + (l-a a )N,i 6 +a a N( — l) 1 ^ = ~ 1 ■ 

Hence, Wgs P ( 2JV )(pi, P2) comes with the coefficient Tr(A) Tr(C) Tr(BID t I). Going in this man- 
ner through all nine cases yields 

(TiiU^AUBU^CUD)) = {^{IU^AUBIUHCUD)) 
= Wg(p 1 ,p 1 )Tr(A)Tr(C)Tr(BD) + Wg(p 1 ,p 2 )Tv(A)Tr(C)Tr(BID t I) 

- Wg(p 1 ,p 3 )Tr(A)Tr( J B)Tr(C)Tr( J D) - Wg{p 2 ,p 1 )Tr{AlC t I)Tr{BD) 

- Wg(p 2 ,p 2 )Tr(AIC t I)Tr(BID t I) + Wg(p 2 , p 3 ) Tr (A* I CI) Tr(B) Tr(D) 

- Wg(p 3 ,p 1 )Tr(AC)Tr(BD) - Wg(p 3 , p 2 ) Tr(AC) Tr(5iT>*I) 
+ Wg(p 3 ,p 3 )Tr(AC)Tr(S)Tr(L>) , 
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and then the claim follows from ([2"T]) . 



A. 4 The general Weingarten integration formulas 

This appendix recollects the general results on the integration over the orthogonal, symplectic 
and unitary group as given in (or can be deduced from) [CSj . |CSt] and |Col] respectively. The 
following notation will be used. For m,n G N denote by J-"(m, n) the set of all maps from 
{1, . . . , m} to {1, . . . , n}. If m is a partition of {1, . . . , m}, then J^m, n) denotes the set of 
those elements of J-"(m, n) which are constant on the blocks of m. Denote by VV(k) the set of 
all pair partitions of {1, ... , k}. In particular, VV(k) — if k is odd. If k — 21 is even, we write 
W(k) 3 m = {{m^(m), n^m)} | v = 1, . . . , /}, where the numbering of the blocks is arbitrary, 
but m v < n v holds for all v. For maps a G F{k, {0, 1}) from {1, . . . , k} to {0, 1}, let us write 
a G A(m, {0, 1}) if for all v holds cr(m„(m)) = 1 — a(n u (m)). 

Now let us begin with the orthogonal case. If N > I, then the space of 0(iV)-invariants in 
V® 21 admits a basis {9 m | m G VV(2l)} given by 

dm = e ^( 1 ) ® ' ' ' ® e «H 2 *) ' 

06^(111,^) 

where as above ej, z = 1, . . . , TV, is the standard orthonormal basis on V = C N which is furnished 
with the symmetric bilinear form b(v, w) = v t w. The Weingarten function Wgo^^m, n) is the 
(m, n)-entry of the inverse of the Gram matrix of the extension of b via b(vi (g) . . .<2>f2j, w i ® • • -® 
W21)) = b(vi,Wi) • . . . • b{v2uW2i) with respect to the basis {9 m \ m G VV(2l)}. For example, on 
VV{2) x VV(2) one obtains Wg 0(7V) = ^, while on VP {A) x PP(4) the orthogonal Weingarten 
function Wg (Ar) is given by the matrix entries 

1 




N(N -1)(N + 2) 

The general orthogonal Weingarten integration formula now reads as follows. 
Proposition 1 Letcp,^ G F(k,N). Then 

I du [ u m,m = ^KroO) MmioCVO w So(jv)(ttv,n) 

JO(N) , t^o/-^ 



i=l m,ngP-p(A:) 



Next let us consider the symplectic case, hence V = C 2N with a as in Section fA. 21 If N > I, 
then the space of SP(2iV)-invariants in V® 21 admits a basis {# m | m G VP (21)}, where 
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with 

v(r),e,j) = { e ;^> ii-™"' ■ 

{ e J/ ( ?7+n (i_ e )) (-1) , it j = n u . 

The Weingarten function Wggp^^^m, n) is the (m, n)-entry of the inverse of the Gram matrix 
of the extension of a (as above) with respect to the basis {9 m | m G VV(2l)}. On VV(2) x W(2) 
one obtains Wg SP ( 2 7v) = ^y, and on VP (A) x PP(4) it is given by (j2"T|) above. 

Proposition 2 Let <f>,ij) G F{k, N), a, /3 G F(k, {0, 1}). Then 



L 



dU ] [ U Hj)+NaU) , m+Nm = Yl W §sp ( 2iv)Kn)(-l)^ a (-W) + «-W' 



SP(2Af) j=1 m ,neVV(k) 



lj-(m,2JV)(0) ^-A(m,2N)(a) 1 f(n,2N) ("0) ^A(n,2N)(P) ■ 



Finally let us turn to the unitary case. Let S^ denote the full symmetric permutation group 
of the set {1, . . . , k}, and e*, . . . , e* N the dual to the basis ei, . . . , ejv on V = C N w.r.t. the 
standard scalar product c. If N > k, then the space of U(iV)-invariants of V® k (g) (y*)® k (with 
the contragredient representation on the dual space) admits a basis {C n \ ti G S&} where 

c n = e ^(^- 1 d)) ® • • • ® ^(t-h*)) ® e h) ® • ■ • ® e 0W • 

4>eT(k,N) 

The Weingarten function Wg U (jv)( (7 ) r ) is the (cr, r)-entry of the inverse of the Gram matrix of 
the extension of the scalar product c w.r.t. the basis {C n \ n G S^}. Then on Si x Si one has 
\Vgu(7v) = -L, while on S 2 x S 2 the unitary Weingarten function is given by 

1 f N -1 

N(N 2 - 1) V -1 N 

Proposition 3 Let <p, ip, <p' , ip' G T{k, N) . Then 

\ dU ] \ U^j)^(j)U^(j)^(j) = Y] Wg u(Ar) {a, t) 6(<f>,a<f>') 6(i>, rij/) . 

Let us conclude with a comment on higher order Weingarten functions. Their definitions 
are given here have the merit that they can be easily motivated in an invariant theoretic 
context. Their calculation, in particular, for higher orders, can be simplified by remarking, 
e.g., that Wg U (_ ZV )((T, r) depends only on the cycle type of ar^ 1 . This leads to sophisticated and 
computationally feasible expansions of Wg at least in the unitary and orthogonal cases |GMj . 
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